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We report a method to locally measure the penetration depth A in a superconductor 
by detecting the diamagnetic response using magnetic force microscopy (MFM). We 
extract A by fitting the height dependence of the levitation force in the Meissner state 
using an analytical model that approximates the MFM tip as a single-domain, trun- 
cated conical shell. We demonstrate on two YBa 2 Cu30 6+a; single crystals with two 
MFM tips that the obtained values agree well with previous results. This approach 
is not affected by the tip width and can be applied to similar but not identical tips. 



PACS numbers: 68.37.Rt, 74.72.-h, 74.25.N- 



The magnetic penetration depth A, determined by the density of the superconducting 
charge carriers, is one of the two fundamental length scales in superconductors.^ Accurate 
determination of A is important for understanding fundamental properties of the super- 
conducting states, such as the order parameter symmetry and the pairing mechanismP^ 
However, the absolute value of A is notoriously difficult to measure. Established methods, 
including muon-spin-rotation,Sl infrared spectroscopy,^ micro-wave cavity techniques^ and 
lower critical field measurements,^ average over bulk samples. Spatially resolved techniques 
may be helpful especially if sample topography, inhomgogeneity, or intrinsic variation are 
suspected. 

Efforts to measure A by magnetic scanning probes include scanning SQUID susceptometryj 
mostly limited by the accuracy in determining the sensor-sample separation z, and magnetic 
force microscopy (MFM) by imaging individual vorticesP^ The quantitative determination 
of A from vortex images is highly non-trivial due to the convolution with the tip structure, 
which requires numerically calculating the convoluted signal based on detailed knowledge 
of the tip magnetic structureP 

In this letter, we report MFM measurements of the absolute value of A by measuring and 
modeling the height dependence of the diamagnetic response in the Meissner state. The 
diamagnetic response is not as badly affected by the finite tip width as lateral imaging, 
allowing us to make approximations that give an analytical description justified for similar 
tips. We demonstrate on two Yf^CusOe+x (YBCO) single crystals that the values we 
obtained with two different tips both agree well with previously published results. 

The measurements were performed in a home-built variable-temperature MFM appara- 
tus. We use high resolution cantilevers (NC-18 from Mikro-Masch), and coat the tip by 
electron beam deposition with an iron film of nominal thickness 40 nm. We measure the 
change of the resonant frequency,^ which is proportional to dF z /dz, where z is normal to 
the cantilever and to the crystal a-b surface, and F is the force that results from integrating 
over the entire tip. The YBCO single crystals were grown by the self-flux method in BaZr03 
crucibles 1 -^ and annealed, with superconducting transition temperature T c f» 56 K, implying 
x ~ 0.56. The samples are platelet shaped, with the face parallel to crystal afe-plane about 
1mm x 0.7mm and thickness 60-80 ycm. The samples were kept for less than two months 
after growth either in a room temperature desiccator or below 77 K ensuring that they were 
fresh. Sample I is fully detwinned, and sample II has twin boundaries separated by a few 
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FIG. 1. Measuring A and A A by MFM in the Meissner state using time-reversed mirror approxi- 
mation, a: z dependence of dF z /dz (grey dots) by MFM at T = 6 K on two YBCO single crystals 
(T c Rj 56 K) and the fit to the truncated cone model (dashed line) from which we extract A as a 
fit parameter. The vertical offset, dF z /dz\ z=00 is and 100 pN///m respectively. Inset illustrates 
the time-reversed mirror approximation: the response of the superconductor can be replaced by an 
image tip reflected over a plane A below the superconducting surface, b: AA(T) = A(T) — X(6K) 
of sample II (grey dots) from T=6 K to 12 K determined independently of the tip model. The 
linear fit (solid line) shows a slope of 0.97 nm/K. 



microns. 

To obtain A, we cool the samples in the absence of magnetic field and measure dF z /dz(z) 
(Fig- (3a)). 

The Meissner levitation force on a magnet from a superconductor occupying half space 
was given in Eq.(2.18) of ref^from which we obtain the force derivative: 
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where r = R + zz, R = RR, Jn is the zeroth order Bessel function of the first kind, and for 
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Approximating G(Xk) by e _2Afc is equivalent to replacing the response of the supercon- 
ductor by an image of the field source mirrored through a plane A below the superconducting 
surface (illustrated in Fig. []Ja)). Under this time- reversed mirror approximation, the change 
in A is nearly equivalent to change in z, allowing us to determine AA(T) = A(T) — X(6K) 
independent of any model of the sensor structure.^^ We obtain from sample II a linear 
AA(T) (Fig. [j](b)) with a slope consistent with previous results,^ as expected from the 
nodal d-wave gap structure of YBCOP 

To extract A, we model the tip as a sharp, single-domain conical shell truncated at the 
distance h . M(r) = M t5(R — a(z + h Q )), where we assume M is along z, a ~ 15° is the cone 
angle, t is the magnetic film thickness, M is the magnetization, and ho is the truncation 
height. We take the approximation of Jo(A;|R — R'|) pa Jq(0) in the limit of a <C 1 and large 
z. From Eq. [T], we obtain: 
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Taking the time-reversed mirror approximation in Eqj2j we obtain: 
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where A = —fi Q (aM tA6) 2 /27r and = 6 2 — 6± = n for our half-coated tips. Since M and 
t may not be known exactly we take A as a fitting parameter. 

We fit touchdown curves in Fig. [IJa) using Eq. |4j fixing ho from scanning electron mi- 
croscopy (SEM) images of the tip (Fig. f2ld-f)) and letting A, X and dF z /dz\ z=OQ vary. In 
the fit, we minimize \ 2 = ZXl — 9F z /dz(z)data/dF z /dz(z) f it ) 2 for z > 0.6/xm (Fig. [2|. We 
obtain A(6iC) = 185 nm and 196 nm in sample I and sample II respectively. Extrapolating 
using AX(T)/T = 0.97 nm/K (Fig. pTb)), we obtain A(0) = 179 and 190 nm, in good agree- 
ment with previous reported values A(0) = 180 ± 20 nm on similar crystals^^feummaries 
in Table [H 

How accurate is the measurement of A? Using the measurement on sample II as an 
example, if we consider only statistical errors, we obtain X a b(6K) = 196 ± 3 nm from 



TABLE I. A of YBa2Cu30g+x {x ~ 0.56) showing good agreement between our results and previ- 
ously published results. The technique, the obtained value of A, the measurement temperature T m , 
and T c are provided. X a b denotes the the average of A along crystal a and b axes (X a b = \AViA&). 



method 


A 


[nm] T m [K] 


T, [Kl 


Ref. 


usR 


A t, 

^ab 


= 175 


1.25 


59 


refP 


lnwpv PT*itir*al ti pi n 


\ >, - 
^ao — 


= 175 ± 6 







ref El 


ESR 


Aa = 


202 ± 22 





56 


refP 


(Gd-doped) 


A fe = 


140 ± 28 








Infrared 


Aa 


= 248 


12 


59 


refP 


spectroscopy 


A 6 


= 183 








MFM 


Kb = 


195 ± 28 


6 


56 




MFM 


A a 6 = 


180 ± 30 


6 


56 





bootstrapping' with 70% confidence interval. 

The systematic errors mainly come from the uncertainty in determining ho and z, and 
the approximation made in the model. The ±20 nm uncertainty on ho leads to ±14 nm 
uncertainty in A. We detect the surface within ±5 nm owing to the abrupt change of 
dF z /dz when van der Waals force dominates over magnetic force (Fig. [ija)). We calibrate 
the scanner using the laser interference pattern from the sample, which gives at most 3% 
error in determining z mostly from the nonlinearity of the scanner. We choose to measure 
fresh surfaces because any non-superconducting layer makes the measured A larger than the 
real value by the thickness of the layer. 

We observe no change in the tip magnetic strength before and after ramping the magnetic 
field to 0.1 T at 30 K, consistent with the tip being mono-domain. The magnetization of the 
tip is presumably aligned along the film. Assuming M along z does not induce systematic 
error in A because the in-plane component gives the same functional dependence of dF z /dz 
as the z component. 

The model ignores the tip width. To estimate the associated error, we release this approx- 
imation by taking the Bessel function in Eq. [l]to the second order: Jq{x) = 1 — 1/Ax 2 + 0(x 4 ) 
and obtain: 

^±( z ) = _g( aMo t) 2 ^fc 4 e- 2fc(A+2) jd6'd9"dz'dz" 
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FIG. 2. a-c: x 2 for the fit in Fig. [ija) for sample II. a: x 2 plotted as the color scale in the 
two-dimensional parameter space of A and A with dF z /dz optimized at each point, b-c: \ 2 as 
a function of A when A is fixed (b) at the value that minimizes x 2 as plotted by the horizontal 
dashed line in (a), and when A is free floating (c) as given by the diagonal dashed line in (a). 
The distribution of x 2 (g re y dots) are from bootstrapping 200 times, d-f: Scanning electron 
microscopy images of the tip used for sample I (d) and for sample II (f), from which we determine 
ho = 120 db 20 nm and 160 db 20 nm respectively. The curvature at the tip apex comes from extra 
deposition of materials due to a sharp edge. 
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Carrying out the integral, the correction on dF z /dz is: 
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The first term on the right hand side of Eq. [5] has the same functional form as the model 
in Eq. |4j so only the second term leads to errors in A, at most 2% owing to the smallness of a. 
Numerical simulations of a truncated cone tip with realistic width quantitatively confirms 
the analytical model of Eq. |4} If the tip is a sharp cone, e.g. ho = 0, the finite tip width 
does not change the functional form of dF z /dz at all. 

Adding all the errors, we obtain A(0) = 190 ± 28 nm for sample II, and for sample I 
A(0) = 179 ± 30 nm. 



In YBCO crystals, A along the crystal a and b axes are not equal, and we measure the 
average A a & = y/X a K- When pinning is weak, we can determine X a /K by resolving the 
positions of individual field-cooled vortices. We demonstrated previously that the Fourier 
transform of vortex positions in a fully doped YBa2Cus07_5 shows an elliptic band, the 
eccentricity of which gives X a / Xb- However, the weak pinning condition may often not be 
satisfiedPI 

To conclude, we demonstrate a method to measure A within 15% error by MFM from 
the height dependence of the diamagnetic response. The essential elements of the method 
are the precise height determination, and modeling the diamagnetic response with a small 
number of parameters. 
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